Spin states in InAs/AlSb/GaSb semiconductor quantum wells 
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We investigate theoretically the spin states in InAs/AlSb/GaSb broken-gap quantum wells by solving the 
Kane model and the Poisson equation self-consistently. The spin states in InAs/AlSb/GaSb quantum wells are 
quite different from those obtained by the single-band Rashba model due to the electron-hole hybridization. The 
Rashba spin-splitting of the lowest conduction subband shows an oscillating behavior. The D'yakonov-Perel' 
spin relaxation time shows several peaks with increasing the Fermi wavevector By inserting an AlSb barrier 
between the InAs and GaSb layers, the hybridization can be greatly reduced. Consequently, the spin orientation, 
the spin splitting, and the D'yakonov-Perel' spin relaxation time can be tuned significantly by changing the 
thickness of the AlSb barrier 

PACS numbers: 71.70.Ej, 72.25.Rb, 73.21.Ac, 73.21.Fg 



I. INTRODUCTION 



very important for many electronic properties. 
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Heterostructures based on InAs, GaSb, and AlSb are one 
of the most promising systems for fundamental physics re- 
search and design of novel devices due to their advantages 
of high electron mobility, naiTow band gap, strong spin-obit 
coupling, and more importantly, the special broken-gap band 
lineup at the InAs/GaSb interface'-^. In the past decades, the 
properties of InAs/GaSb broken-gap supperlattices and quan- 
tum wells (QWs) were investigated both theocratically and 
experimentally ^''-^'^'^'^-^-^' ' 'i '2, 13. 14 ^ These studies show that 
due to the overlap of InAs conduction band and GaSb valence 
band, the energy dispersion may exhibits an anticrossing be- 
havior at a finite in-plane wave vector 9^ q3.4.5,6,7,8,9,io ^ rpj^g 
tunneling between InAs and GaSb layers opens a mini hy- 
bridization gap, which was observed experimentall y ' ' ^i ' 3, i4 ^ 
In addition, the electrons in GaSb can move across InAs/GaSb 
interface into InAs layer, forming a two-dimensional electron 
gas in InAs side and a two-dimensional hole gas in GaSb side, 
which is promising for observing the Bose-Einstein conden- 
sation of excitonsiiii^ii^. In practical applications, there have 
been many proposals for electronic and optical devices utiliz- 
ing the unique characteristics of InAs/AlSb/GaSb system such 
as resonant tunneling structure a'^i'^ , infrared detectorsi^, and 
interband cascade laser diodes^". Recently, the spin-related 
properties of InAs/AlSb/GaSb system also attracted much in- 
terest. For instance, there have been a number of spintronic 
device proposals, including the Rashba spin filter^', the spin 
field effect transistor^^, and the high-frequency optical mod- 
ulator utilizing the spin precession^^. Moreover, InAs/GaSb 
QW, like HgTe/HgCdTe QW, is another possible candidate to 
demonstrate the intrinsic spin Hall effect^^ and quantum spin 
Hall phase^^ due to the inverted band structure. 

In the previous works, the electron-hole hybridization in 
InAs/GaSb broken-gap QWs have been well studied2ii2i2^. 
It was shown that the "spin-up" and the "spin-down" states 
are affected differently by the hybridization'' "'. This implies 
that the spin states near the hybridization gap (where strong 
hybridization occurs) in InAs/GaSb QWs should be quite dif- 
ferent from those in conventional semiconductor QWs. Since 
the Fermi level of the undoped InAs/GaSb QW lies inside the 
hybridization gap^^, the unusual spin states nearby would be 



In this paper, we investigate theoretically the spin orienta- 
tion (i.e., the expectation value of the Pauli operator <t = 2S 
for the electron spin S in an eigenstate), the zero-field spin 
splitting, and the D'yakonov-Perel' (DP) spin relaxation— 
in InAs/AlSb/GaSb QWs based on the Kane model. The 
charge transfer induced internal electric field is taken into ac- 
count by solving the Kane model and the Poisson equation 
self-consistently. For InAs/AlSb/GaSb QWs, the spin orien- 
tation (or spin states) are significantly modified by the hy- 
bridization between the conduction band and valence band, 
which can be tuned by changing the thickness of AlSb bar- 
rier. The spin orientation can be measured indirectly from 
the Faraday rotation, and indeed offers us a physical picture 
and can lead to some novel effects such as the persistent Spin 
Helix2S. Using the grating technique and Faraday rotation, 
the persistent rotation of the spin orientation, i.e., spin he- 
lix of electron was observed in GaAs QW^". The Rashba 
spin-splitting (RSS) exhibits oscillating features as a function 
of in-plane wave vector in InAs/AlSb/GaSb QW— near the 
hybridization gap. The spin relaxation time, which are ob- 
tained from the perturbation theory^^, shows that the unusual 
spin-splitting could lead to several peaks with increasing the 
Fermi wave vector. These unusual features all come from 
the strong electron-hole hybridization beyond the single band 
model with linear Rashba spin-orbit interaction^^. Interest- 
ingly, the hybridization is very sensitive to the thickness of the 
AlSb barrier, evidenced by a rapidly decreasing hybridization 
gap with increasing AlSb barrier thickness. Consequently, all 
the spin-related properties, including the spin orientation, the 
spin splitting, and the DP spin relaxation time, can be tuned 
by varing the thickness of AlSb barrier. The property of tun- 
able spin states in InAs/AlSb/GaSb QWs might be useful in 
designing new spintronic devices. 

This paper is organized as follows. In Sec. |II] we describe 
the theoretical method based on a self-consistent calculation 
combining the Kane model and the Poisson equation. In Sec. 
HID we present the numerical results for the band structure, the 
spin orientation, the spin-splitting and DP spin relaxation time 
in InAs/AlSb/GaSb QWs. In Sec. IV we give the conclusion. 
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II. THEORY 

We consider an InAs/AlSb/GaSb broken-gap QW grown 
along the [001] direction [see Fig. [I](a)]. We choose the axis 
X, y, and z to be along [100], [010], and [001], respectively. 
Within the envelope function approximation, the Kane model 
is a good starting point for systems with strong electron-hole 
hybridization like InAs/AlSb/GaSb QWs. As InAs conduc- 
tion band overlaps with GaSb valence band, electrons could 
transfer from GaSb layer into InAs layer. The charge redistri- 
bution induces an internal electric field, which can be evalu- 
ated from the Poisson equation. Generally, one needs to know 
the charge density distribution to solve the Poisson equation, 
and the charge density distribution is in turn determined by the 
electron wave function. Therefore a self-consistent procedure 
is needed to take the charge-transfer effect into account in this 
system-''^. 

In Sec. Ill Al we discuss the Kane model and the self- 
consistent calculation method. To consider the spin states in 
broken-gap QWs, we shall give an explicit definition of "spin- 
up" and "spin-down" states by classifying the eigenstates of 
the Kane model. This is discussed in Sec. IIIBI 



A. Hamiltonian and calculation method 



can be generalized to describe heterostructures by ordering 
the momentum operators with respect to material parameters. 
By choosing the following set of basis functions. 
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Following the new envelope function theory^^, the Kane 
model which describes the bulk zinc -blende semiconductors 



the Kane Hamiltonian for zinc-blende crystals near the F 
point is 
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where In Eq. (O, k - {k\\^,-idld.), k± = + iky, and {kaykp} = 

(kaykn + knykaj 12 (a,p - x,y,z)- Here the in-plane momen- 

^ ' ^ ' tum as a constant of motion has been replaced by its eigen- 

p _ j7 ^ ]^ ]^ (2h) v^lus ^11 ^iid we have neglected the bulk inversion asymmetry 

~ ^ 2mo (DresseUiaus effect) since it is small compared to the struc- 

^2 ture inversion asymmetry (Rashba effect). The strain effect is 

Q - - — (kxjikx + ky-yiky - 2k{y2kz), (2c) also neglected because the lattice mismatch between InAs and 

" GaSb is less than 1%. 

L = /- [k^y^k.], (2d) 

M = - [k,y2kx - k,y2ky - 2/{fe,y3^v)] , (2e) 

ZniQ 

U = -^Pok„ (2f) 
V3 

y = -^Pok-. (2g) 
V6 
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TABLE I: The Kane parameters used in our calculation. These pa- 
rameters are calculated from the Luttinger parameters obtained from 
Ref.!40 (see the formulism below Table|IJ. The value of Kane energy 
Ep = linoPl/h^ is taken equal to 22.5 eV for each layer material. 





InAs 


GaSb 


AlSb 


E, (eV) 


0.417 


0.812 


2.386 


£v (eV) 


-0.417 


0.143 


-0.237 


A(eV) 


0.39 


0.76 


0.676 


A, (eV-nm^) 


-0.26 


0.09 


-0.06 


n 


2.01 


4.16 


2.04 


72 


-0.49 


8.18 


-0.38 


73 


0.21 


1.38 


0.40 


e 


14.55 


15.69 


14.4 



"The relation between Kane parameters and Luttinger parameters are yi = 
'I lEg ' 72 - ?2 6£g ' 71 - 73 6£j ■ 
''Ac is correlated to electron effective mass m* by A^. = ^-r - + 



The band parameters are assumed to be a piecewise func- 
tion along the growth direction, 

r fe) = X ^' [® " ^'^ " ® " ' 

where 0(z) is the Heaviside step function, z, is the i-th in- 
terface of this system, and y' is the bulk band parameters 
of the i-th layer. These band parameters, including Eg, Ey, 
A, Ac, 71, 72, and 73, can be derived from the Luttinger pa- 
rameters and the electron effective mass"*". They are given in 
table H] The valence band offset of two neighboring materials 

is A^;, = £■;, -£■;+'. 

For thick enough InAs and GaSb layers, the lowest con- 
duction subband in InAs layer overlaps with the highest va- 
lence subband in GaSb layer. Electrons could transfer from 
GaSb to InAs, inducing an internal electrostatic potential 
Vin (z). Therefore, the total Hamiltonian becomes H{k\\) = 
Hk{k\\)-eVin (z). The subband dispersions and the correspond- 
ing eigenstates are obtained from the Schrodinger equation 

7/(*ll)|^,(^ll)> = £,(^11) 1^,(^11)), (4) 

where s is the index of the subband and |*I'5(^||)) - 
exp(/^ll'P)[i^j(z), ^fjCz), fl{z)Y is the envelope function. To 
solve the Schrodinger equation, we expand by a series of 
plane waves, 

^P,'(z) = — ^ c^,„ exp(il„,z), (5) 

m=-N 

where km - Imn/L and L is the total length of the structure 
(in this work, L = IL^f^,^ + Li„as + L'^lf,'' + LcaSb. and U^^^ = 
2Lj„As)- By choosing a moderate A^, one can also avoid the 
spurious solutions that may occur in the Kane model. In our 
calculation, A^ » 25 is good enough to get convergent results. 

The internal electrostatic potential V,,, (z) is determined by 
the Poisson equation 



-s(z)-Vi„iz) = - Ipciz) + Phiz)] , (6) 
az dz 

where Pe{z) and p/,(z) are, respectively, the charge density due 
to electrons and holes and e(z) is the static dielectric constant. 
Peiz) and Phiz) can be derived from the envelope functions: 

P^(^) = -TTTi Z r Z V'nizf fFiEs)dh, (7) 

Phiz)^-:^y{ V \ip:iz)\\l-fFiEMdh, (8) 

where fpiEs) is the Fermi distribution function. The summa- 
tions 2«=i.2 and 2„=3,4,....8 runs over the electron components 
01, 02 and hole components 03, ■ ■ ■ ,08 respectively. The sum- 
mation Ej includes all the subbands which show anticross- 
ing behavior. For simplicity, we take T = K and the axial 
approximation22i2^ in the self-consistent procedure. Since we 
want to give a clear picture about the effect of the hybridiza- 
tion on the spin states, we only consider the anticrossing oc- 
curs when the lowest conduction subband meets the highest 
valence subband, which limits to the two cases Li„as ^ 14 nm 
at LcaSb - 10 nm and LcaSb < 14 nm at Li„as - 10 nm. 

The Fermi level Ef k determined by the charge neutrality 
condition. 

r \peiz)+Phiz)]dz^Q. (9) 
Jo 

The Fermi level Ef obtained by Eqs. (jTji-® locates closely 
above the mini gap^^-^'*. Combining Eqs. (I4ll-(|9]l, we can 
do a self-consistent iteration that eventually yields the inter- 
nal electrostatic potential y,„ (z). Once we have V,,, (z), the 
subband dispersions and electronic states can be obtained by 
solving Eq. (|4]i. 

B. The classification of the spin states 

In the single-band Rashba model for the electron (i.e., 0i 
and 02 components only), there is a well-defined spin quanti- 
zation axis By, - e^y- e, perpendicular to both the wave vector 
direction e^^^ and the QW growth direction e,. The spin orien- 
tation of any eigenstate always has a magnitude of unity and 
is either parallel (called "spin-up" eigenstates) or anti-parallel 
(called "spin-down" eigenstates) to Cjr^. For the 7 = 3/2 hole 
system (i.e., components 03, ■ ■ ■ , 06) or a hybridized electron- 
hole system/-', however, such a quantization axis does not ex- 
ist. As we shall show in Sec. |llll for a given conduction 
subband, the spin orientation may change its magnitude and 
direction (up to 180 degrees) with increasing k\\. The absence 
of a well-defined spin quantization axis make it impossible 
to classify the "spin-up" and "spin-down" states by their spin 
orientations relative to this quantization axis. 

In the following, we divide the Kane Hamiltonian Hkik\\) 
in Eq. (O (note that bulk inversion asymmetry has been 
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neglected) as the sum of a dominanting part Hax{k\\) with 
axial symmetry and a small cubic part Hcutikw) with cubic 
symmetry' 



and 
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Hk(k\\) - Hax{k\\) + Hcub(k\\). 



(10) 



The axial part Hax(k\\) obtained from Eq. ([TI]i by replacing 
72 and 73 with 7 in the term M and the cubic part Hcub(k\\) = 
Hkik\\)-Hax{k\\) is the difference between Hk{k\\) and Hax{k\\). 

The axial part Hax{k\\) can be transformed into a block- 
diagonal form by choosing a new basis set. Similar trans- 
formation has been reported in dealing with the four-band22i 
and six-band models^, but the cubic part was neglected in 
these works. In this paper the transformation is extended to 
the eight-band model and the cubic part is also included. Let 
^11 - fc||(cos (fi, sin if) and the new basis set is (/)-dependent, 
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Under this basis, Hax(k\\) 'J^ax(k\\) is block-diagonalized 
and Hcub{k\\) — > 'Hcub(k\\) has a simple dependence on the az- 
imuth ifi. 



'Hax(h) : 



'J-(cub(k\\) = 



•H-(^ll) 
-^+(^11). 

y[(A:||)cos49J -/y[(A;||) sin4vp'* 
/^(A:||)sin4v3 -^(A;||) cos4<^ 



(12) 
(13) 



where 



VJy - V2f/' + y f/' + V2y' 



V3y' -p-Q 

y/2U' + V -U + M' 

-u' + V2y' 



L' + M' 

-p + e 

-2gT VIl' 
V2 



L'±2M' 
V2 

-2Q± V3L' 
V2 

-P- A 
(14) 



2ot() 



-A7 





1 V2 

10 

V2 



(15) 



L' = i— fc||(73fc + k,y3), 



M' 



2mo 
2mo " 



U' = — Po^z, 
V3 

1 

A7= ^(72-73). 



(16a) 

(16b) 
(16c) 
(16d) 

(16e) 
(16f) 



Using the new basis set, we can give an explicit classifi- 
cation of all the eigenstates into "spin-up" and "spin-down" 
states, similar to the single-band model. 

In the absence of the cubic term ^ciib{k\\) (so-called axial 



approximation), the eigenstates \^°'l^{k\\)) and eigen-energies 



^s°±'(^ll) (for the i-th subband) of the total Hamiltonian 
H{k\\} - 'Haxik\\) - eVin (z) are determined by 

[9ijk\0 - eVin(z)] \'i't::!(k\0) = £?±'(^||) [i'ifikii)) , (17) 
with the Rashba spin splitting in the .s-th subband 



^{ax) 



■{ax). 



(18) 



Obviously, the eigenstates are automatically classified into 



two classes: the "spin-down" states \^"^{k\\}) consisting of 
the components |5(-)>, \HH(-)), \LHi-)), and \S 0(-)) and 
the "spin-up" states |*I'j'';^'(^||)) consisting of the components 
\S(+)), \HH{+)), \LH(+)), and \S0(+)). Then one can define 
the spin orientation of an arbitrary eigenstate \^^s.±(k\i)) as^^^ 

0"Um,,^ = (*pS^(^I|)| K±(k\\)) , (19) 

where L(^) = [l,xi(p),l.yitp),1.7i(p)] are 8 x 8 matrices (see 
Appendix lAb for the Pauli operators cr = 25 in the new basis 
(El). 

In the presence of the cubic part 'Hcubik\\), the off-block- 
diagonal terms +i^sm4(f will in general couple the "spin- 
up" states with "spin-down" states, unless points along a 
high symmetry axis satisfying sin4(/? = 0. As a result, the 
exact eigenstate |*I's(^||)) of the Kane Hamiltonian is in gen- 
eral a mixture of "spin-up" and "spin-down" states. For con- 
venience, however, we still classify the exact eigenstates as 
"spin-up" or "spin-down" states according to the dominant 
component. 

Now we discuss the electronic structure of InAs/GaSb and 
InAs/AlSb/GaSb QWs in the axial approximation ['Hcubik\\) = 
0] and the modification due to the cubic correction 'Hcub{k\\). 
In the axial approximation, the electronic structure shows 
three distinct features: 
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Electron-hole anticrossing. If the InAs and GaSb layers 
are thick enough, the lowest conduction subband (CBl) 
anticrosses with the highest valence subband (VBl) at 
a critical in-plane wave vector k^, evidenced by the 
switch of the dominant component (a) from \HH(+)} 

to \S(±)) for the CBl subband state I^Pcbi ±('^ll)) 
(b) from |5(+)> to \HH(+)} for the VBl subband state 



(ax) 
VBl,± 



when is increased across ka- 



2. The spin orientation of is always along e^. - 

X i.e., in the QW plane and perpendicular to 
the in-plane wave vector k\\, in agreement with the pre- 
scription of the single-band Rashba modeb^^ "*^ (this can 
be verified using Eqs. ( IA8al i-( |A9cl ) in Appendix [A]). 
But the projection of the spin orientation varies with 
and may even change sign, which is different from the 
single-band Rashba model. 

3. The Rashba spin-splitting of the CBl subband exhibits 
an oscillating behavior [containing zero spin-splitting 
(spin degeneracy) points] due to the electron-hole anti- 
crossing. 

In the presence of the cubic correction l-icubikw), the follow- 
ing additional features are introduced: 

1 . The band structure and spin-splitting as a function of k\\ 
display a C4y symmetry. 

2. The spin orientation of the exact eigenstate |*I's(^||)) de- 
viates from = x e, (although it still lies in the 
QW plane) due to hybridization of spin-up and spin- 
down states. The maximum deviation arises when the 
azimuth angle (/? of k\\ satisfies sin(4^) = +1 [such that 
the spin mixing matrix element in Eq. (fTST i reaches the 
maximum] and at one of the spin degeneracy points. 

In Sec. Uni all these special properties in InAs/GaSb and 
InAs/AlSb/GaSb broken-gap QWs are demonstrated by our 
numerical results. 
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FIG. 1 : (Colour online) (a) Band profile and the probability density 
distribution of CB 1(±) (blue line), VB 1(±) (red line) states in a 10-2- 
10 nm InAs/AlSb/GaSb quantum well at fcy = 0; (b) the amplification 
of the dashed square area in panel (a). The black solid, blue dashed, 
red dashdot and green dotted line in panel (b) denotes the results of 
I^Aisb = nm, 1 nm, 2 nm and 5 nm respectively. 
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III. NUMERICAL RESULTS AND DISCUSSION 

In Fig. [T] we show the self-consistent band profile of an 
InAs/AlSb/GaSb broken-gap QWs. In this QW, the concen- 
tration of electrons transferred from GaSb layer to InAs layer 
is found to be of the order 1.61 x 10" cm"^. The resulting 
electronstatic potential Vi„{z) induces a 8~14 meV downward 
(upward) bending of the InAs conduction band (GaSb valence 
band) near the interface and causes a slight shift of the sub- 
bands (see Fig. [I](b) ). Note that we denote the upper two 
branches of anticrossed subbands as CBl{+) and the lower 
two as VBl{+) (see Fig. |2]i, which are different from the pre- 
vious works^'"'. 

In Fig. |2] we plot the band structure of InAs/AlSb/GaSb 
QWs with different thicknesses of middle AlSb barrier (The 
thicknesses of InAs and GaSb layers are fixed at 10 nm). 
Due to electron-hole hybridization, the CBl and VBl sub- 
bands exhibit a strong anticrossing behavior and open a mini 



FIG. 2: (Colour online) The subband dispersions of InAs/AlSb/GaSb 
quantum well structures with dilferent thicknesses of AlSb layers: (a) 
LAisb = (b) Laisii = 1 nm (c) Labs = 2 nm (d) LAisb = 5 nm. The 
thicknesses of InAs and GaSb layers are fixed at 10 nm. The Fermi 
energy denoted by the dashed lines lies at the bottom of the lowest 
conduction band, in agreement with previous works----. In panel (a) 
and (c), the cross points between CBl (VBl) and the purple (green) 
dashdotted lines mark the constant energy contours drawn in Fig. [3] 

(Fig. a. 



hybridization gap at a finite k\\, consistent with the previous 
works^i^ii2i2ii2. The resulting strongly-hybridized states near 
the gap may have significant contributions to the spin-related 
properties of the broken-gap QWs, since the Fermi energy 
(the dashed lines in Fig. |2| locates nearby22i21. From Fig. 
|2](b)-(d), we can see that by increasing the thickness of the 
AlSb barrier layer, the anticrossing between CBl and VBl 
is gradually weakened and the hybridized gap is noticeably 
narrowed due to the tunneling between CBl and VBl is sup- 
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FIG. 3: (Colour online) The spin orientations on the constant energy 
contours of CS1(+) subband of (a) a 10 nm InAs/10 nm GaSb QW; 
(b) a 10 nm InAs/2 nm AlSb/10 nm GaSb QW. Note the dominant 
components of the states of the insidest contours are \HH(+)), while 
that of the other contours are |S (+)). 



FIG. 4: (Colour Online) The spin orientations on the constant energy 
contours of VBl subband of (a) a 10 nm InAs/10 nm GaSb QW; (b) 
a 10 nm InAs/2 nm AlSb/10 nm GaSb QW. The blue and red lines 
represent VBl(+) and VBl(-) branches respectively. 



pressed significantly. Meanwhile, the spin splitting of each 
subband is greatly reduced due to the decreasing structural 
inversion asymmetry in the InAs/AlSb/GaSb QW compared 
with the InAs/GaSb QW. 

In Fig. [3] (a) we show the spin orientations of the eigen- 
states on different constant energy contours of CB1{+) sub- 
band in InAs/GaSb QW. The spin orientations on the con- 
tour of CBl(-) subband are antiparallel to those of CBl{+) 
subband and are omitted in the figure for brevity. In the k- 



linear Rashba model, the spin orientations of an eigenstate 
is along Cz = Sk^^ x e^, i.e., along the tangent direction of 
the circular energy contour in the QW planed. However, 
this property no longer holds for InAs/AlSb/GaSb broken- 
gap QW. The spin orientations deviate strongly from the tan- 
gent direction Ck,, x C;, unless points along high symme- 
try directions (such as (100) and (110)) satisfying sin 4^ = 
0. This comes from the hybridization between the "spin- 
up" and "spin-down" states, as discussed in the previous sec- 
tion. Therefore we can see the spin orientations on the contour 
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FIG. 5: (Colour online) The components of states in CB1(±) subband 
as a function of k^^ along tp = n/S (solid line) and = n/4 ([110]) 
direction (dashed line) in a 10 nm InAs/10 nm GaSb QW (a), (b) and 
a 10 nm InAs/2 nm AlSb/10 nm GaSb QW (c), (d). The red arrows 
indicate the maximum hybridization point. 



FIG. 6: (Colour online) The projection of spin expectation values of 
|l'csi.+ ) and l^vsii) states along the orthogonal direction of k\\ in 
InAs/AlSb/GaSb quantum well structures with different thicknesses 
of AlSb layers: (a) LAisb = (b) Lab* = 1 nm (c) Laks = 2 nm (d) 
I^Msb = 5 nm. The thickness of InAs and GaSb layers are fixed at 10 
nm. 



£■ = 130 meV deviate the most heavily because this contour 
is nearest to the maximum hybridization point, and when 
lies in the directions if = 7r/8, 37r/8 (sin 4^ = +1), the max- 
imum hybridization occurs. When we insert an AlSb barrier 
between InAs and GaSb layers, the hybridization is strongly 
reduced [Fig. Ob)]. Thus the deviation of the spin orientation 
from the tangent direction is very small. The results of Fig. [3] 
imply that one can tune the spin orientations near the Fermi 
level by changing the thickness of AlSb barrier in the middle 
of InAs and GaSb layers. 

Fig. |4] (a) exhibits the spin orientations for the states 
on different constant energy contours of VBl subband in a 
InAs/GaSb QW. The energy contours of VBl subband show 
a very complicated behavior and a strong anisotropy in [100] 
and [110] directions due to the complicated band structures 
(see Fig. |2]i. From the figure, one can easily find a C4,, group 
symmetry, which comes from cubic symmetry of the crystal. 
For £■ = 113 meV, we can find two pairs of contours. The 
states are electron-like for the inner but hole-like for the outer 
pair. Due to the large spin splitting, the shape of energy con- 
tour for "spin-up" and "spin-down" state looks very different. 
If we insert an AlSb barrier between InAs and GaSb layers, 
the spin-splitting between VBl{+) and VB1{-) is greatly re- 
duced (see Fig. (HJi (b)). 

In order to demonstrate the hybridization of the "spin- 
up" and "spin-down" state, we plot the components |5(+)), 
\HH(±)), \LH(±)) of states in CB1(+) subbands in Fig. |5] 
From this figure, one can see that the components \HH(+)) 
and |5(+)) varies significantly when increase ^:|| over the an- 
ticrossing point ka"^'^^'^^. Interestingly, at ip - 7r/8 direction, 
a strong hybridization between the "spin-up" state and "spin- 
down" state in the CBl subbands occurs. This feature can be 
proven by the peak of "spin-down" component and the dip of 



"spin-up" component at k\i - 0.145 nm"' in the CBl(+) sub- 
band. A similar behavior appears in the CBl{-)) subband. 
For along [110] {sm{4(f) = 0), there is no hybridization in 
the CBl subbands and leading to pure "spin-up" and "spin- 
down" states . By inserting an AlSb barrier between InAs 
and GaSb layers, the hole-like to electron-like transition of 
|*I'cBi,±) states become more emergent, and the hybridization 
between "spin-up" or "spin-down" state is reduced. 

As a result of the dominant component transition of the 
|*iI'cBi,±) and I^FvBi^-t) states when k\\ sweeps across the an- 
ticrossing point ka, the spin expectation value magnitude (S) 
change correspondingly. In Fig. (jSjl we display the change 
of (E) for |*I'cBi,±) and |*I'vbi,±) states as a function of ^y. 
(Z) can be defined by projecting the vector (L) onto e^, with 
^2 = X ^; is the unit vector of the in-plane direction per- 
pendicular to ^||. A sudden change of (S) appears in Fig. ^ 
when sweeps across the anticrosssing point so that the 
main characteristic of |*I'cbi,±) (|*1'i/bi,±)) states change from 
hole-like (electron-like) to electron-like (hole-like). In addi- 
tion, we find sign reversals occur for (S) near the anticross- 
ing points, which means the spin orientations do not maintain 
the same direction. This leads to the failure of recognizing 
the "spin-up" and "spin-down" branches simply by their spin 
orientations. Therefore we should classify the different spin 
states in a new set of basis functions as discussed in Sec. IIIBI 
Increasing the thickness of middle AlSb barrier, i.e., weaken- 
ing of the interlayer coupling between InAs and GaSb layers, 
makes the smooth variation of (S) more and more sharp. 

Besides the spin orientations in InAs/AlSb/GaSb QWs, 
it is interesting to discuss the zero-field spin-splitting in 
these QWs because it can be directly measured from the 
experiments"*^. Therefore, we plot the Rashba spin-splitting 



8 



. (a) 


Lj,5, = 0nm . 




VB1 




\ /cBl\ 


AE = V H 



(c) 



: 2 nm 




(b) L 


«sb = 1 nm 








(d) L 


Aisb = 5 nm 






: 





> 

2 § 
UJ 

< 



1 

> 

E_ 

HI 

< 




0.6 0.4 0.2 0.2 0.4 0.6 0.6 0.4 0.2 0.2 0.4 0.6 

[110] k,,(nm-') [100] [110] k,, (nm ') [100] 





b) 


y-ust, = 1 nm 


:r 


V 





A (c) 






kisb = 2 nm 


■ 11/ 









LA,sb = 5nm 





0.1 0.2 0.3 0.4 0.5 0.6 0.1 0.2 0.3 0.4 0.5 0.6 

kp (nm"^) kp {nm"^) 



FIG. 7: (Colour online) Rashba spin-splitting of CBl (blue line) and 
VBl (red line) subbands in InAs/AlSb/GaSb quantum well structures 
with different thiclmesses of AlSb layers: (a) L^/si = (b) Laisi, = 1 
nm (c) L^isb = 2 nm (d) Laisii = 5 nm. The thickness of InAs and 
GaSb layers are fixed at 10 nm. The red arrows mark the cross point 
between the curve of RSS and the dashdotted line AE = 0. 



FIG. 8: (Colour online) Dyakonov-Perel spin relaxation time of CSl 
subband as a function of Fermi wave vector kf in InAs/AlSb/GaSb 
QWs with difi"erent thiclmesses of AlSb layers: (a) Lab* = (b) 
i-A/i* = 1 nm (c) LAisb = 2 nm (d) L^isb = 5 nm. The red arrows mark 
the resonant peeks corresponding to the spin-degeneracy points. 



(RSS) of CBl and VBl subbands as a function of the in-plane 
momentum in Fig. Q (a). From the figure, one can see a 
valley and sign-reversal occurs in the RSS of CBl subband, 
leading to the oscillating behavior. This anomalous behavior 
arises from the difference between the anticrossing point be- 
tween the CBl{+) and VBl(+) subbands and that between the 
CBl(-) and VBl{-) subbands. The decrease of RSS appear- 
ing at large ^:|| is caused by the weakening of the conduction- 
valence band-coupling for carriers with large momentum, i.e., 
large kinetic energy or large effective bandgap"*^. In Fig. |7] 
we have marked the AE = (spin degeneracy) points with 
red arrows. By comparing to Fig. |5] we find these points 
actually lead to the maximum hybridized points in Fig. |5] 
The splitting of VBl subband is much larger than that of CBl 
subband. This reflects the fact that the spin-orbit coupling in 
valence band is much stronger than that of conduction band. 
An extremum appears in the RSS of VBl subband near the 
anticrossing point. Fig. |7](b)-(d) shows the RSS of QWs with 
a AlSb barrier inserted between InAs and GaSb layers. When 
introducing an AlSb barrier into InAs/GaSb QW, the asymme- 
try at the left and right interfaces for InAs and GaSb layers is 
compensated, so the RSS of InAs/AlSb/GaSb QWs decreases 
greatly with increasing the thickness of AlSb layer The val- 
ley in RSS of CBl subband becomes sharper because the an- 
ticrossing behavior between the CBl and VBl subbands is 
heavily weakened as the AlSb barrier thickness increases. In- 
terestingly, the anticrossing behavior seems like a crossing for 
the thick middle AlSb barrier, e.g., Laisb = 5nm. 

Since the spin-splitting and spin states in InAs/GaSb and 
InAs/AlSb/GaSb broken-gap QWs are very different from that 
in conventional semiconductor QWs, the spin-related proper- 
ties in these QWs should manifest a distinct feature as a con- 



sequence. As an example, the DP spin relaxation time of CBl 
subband in InAs/GaSb and InAs/AlSb/GaSb broken-gap QWs 
is calculated by taking axial approximation and based on the 
perturbation theory^. This theory demonstrates that the DP 
spin relaxation rate tZ^ oc oc (AEJi^j)', where O is the 
spin-obit coupling induced in-plane effective magnetic field 
and proportional to the spin splitting A/s^^j . The perturbation 
theory gives a clear physical picture about the DP spin relax- 
ation, that the DP spin relaxation time would show resonant 
peaks when the spin splitting vanishes. As shown in Fig. |8] for 
a momentum relaxation time Tj, =0.1 ps, we find the DP spin 
relaxation time in these QWs varies from 1 ps to 10^ ps with 
different kf, and the DP spin relaxation time t- exhibits an ob- 
vious oscillating behavior The resonant peaks (marked with 
the red arrows), actually corresponds to the AE = in (spin 
degeneracy) points Fig. Q From panel (a)-(d), we can see the 
spin relaxation time in InAs/AlSb/GaSb QWs is very sensi- 
tive to the thicknesses of middle AlSb layer The oscillating 
and large-scale variation features of DP spin relaxation time 
in InAs/GaSb/ AlSb QWs are dramatically different from that 
in conventional semiconductor QWs. We suppose the unique 
features of spin relaxation time InAs/GaSb/ AlSb broken-gap 
QWs could provide us an interesting way to manipulate the 
evolution of electron spins . 

In Fig. |9]we show the mini gap as a function of the thick- 
ness of the middle AlSb barrier in InAs/AlSb/GaSb QWs. The 
mini gap describes the degree of electron-hole hybridization. 
Fix the thickness of InAs and GaSb layer of InAs/AlSb/GaSb 
QW and increase the thickness of middle AlSb barrier, the 
mini gap decrease rapidly, this feature was akeady demon- 
strated experimentally in Ref. [l^ in which a 1.75 meV 
mini gap has been measured for a 15 nm-10 nm InAs/GaSb 
QW and a 1.75 meV mini gap for a 15 nm-1.5 nm-15 nm 
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FIG. 9: Mini gap as a function of the thickness of AlSb barrier in 
InAs/AlSb/GaSb QWs with (a) LgaSb fixed at 10 nm, (b) L,„as fixed 
at 10 nm. 



between the InAs and GaSb layers, the degree of hybridiza- 
tion can be tuned heavily, the mini hybridization gap and spin 
orientations on the Fermi surface can be changed as a conse- 
quence. Our theoretical calculation is interesting both from 
the basic physics and potential application of the spintronic 
devices based on this novel broken-gap semiconductor QW 
system. 



Acknowledgments 

This work was supported by the NSFC Grant Nos. 
60525405 and 10874175, , and bilateral program between 
Sweden and China. 



APPENDIX A: EIGHT-BAND SPIN MATRICES 



InAs/AlSb/GaSb QW. The mini gap measured from the ex- 
periments in Refs, 1 llll2lll3l are 4, 7 and 2 meV respectively. 
Based on our calculation, the mini gap is found to be 0~4 
meV, the order agree with the experiments. In general, we 
can see the mini gaps reduce to zero when the thickness of 
the AlSb barrier become larger than 5 nm. This is because 
the tunneling between InAs conduction band and GaSb va- 
lence band is greatly suppressed by the AlSb barrier, so the 
electron-hole hybridization is restricted. In addition, if we fix 
AlSb barrier, and increase the thickness of InAs layer or AlSb 
layer, the mini gap decrease too. This is because the confining 
energy is reduced as the thickness of InAs layer or AlSb layers 
increases, and the anticrossing point is then moved towards a 
higher k\\, which has less conduction-valence interband cou- 
pling strength and forms a smaller mini gap. As the mini gap 
and the hybridization degree change, the spin-related prop- 
erties, including the spin orientations, spin splitting, and DP 
spin relaxation time change consequently. Therefore we cer- 
tainly find a method to tune the spin states in InAs/AlSb/GaSb 
QWs, which might be taken advantages in designing spin- 
tronic devices. 



IV. CONCLUSION 

We have investigated theoretically the spin orientation, 
spin-splitting, spin relaxation in InAs/AlSb/GaSb broken-gap 
QWs. We found the spin states in these broken-gap QWs are 
very different from that in conventional semiconductor QWs. 
The spin orientations deviate away from the tangent direction 
of the energy constant surface and the RSS of the anticrossed 
CBl subband in InAs/AlSb/GaSb QWs exhibits a nonlinear 
and oscillating behavior The deviation of spin orientation 
comes from the strong hybridization between different spin 
states, and the oscillating behavior of RSS is a result of the an- 
ticrossing of energy dispersions according. The distinct prop- 
erties of RSS in InAs/AlSb/GaSb QWs lead to an oscillating 
behavior of DP spin relaxation time obtained from the per- 
turbation theory. By changing the thickness of AlSb barrier 



Following the derivation of R. Winkler—, we can obtain the 
form of eight-band spin matrices vector in the basis set (fThl i 
by E' = cr^horh, where cr - (cr^, cr^, crS) is the vector of Pauli 
spin matrices, and h„rb refers to the orbital part of the set basis 
function (fThl i. The components S^, S' of E' can be written 
as: 



Z' = 
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Note that the set of basis functions we used are a little different 
from the basis used in Ref. 42, so the form of spin matrices 
(IAlb -( IA3b are different from these in Ref. Using Eq. (fTTT) . 
we can transform E'^., EJ,, I.', into Sv(^), Ey((/;), T-zif), which are 
the components of eight-band spin matrices in the new basis 
set. S v(^), Ey((/;), can be written as 



I.y{(p) = 



S sin (fi iC cos ip 
-iC^ cos ip Dsmip 

-S cos (fi iC sin 
-iC^ sin (fi -Dcos(f 



[r 

where S, C, £), ;F are 4x4 matrices 
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For a pure "spin-down" state |*P5,_(^||)), the expectation value 
of 'Ex{(f), 'Ly{(p), E;(v?) can be evaluated by 



<2,(^ll)>,,_ =<»l',,_(fc||)|2:,(^)|T,,_(^||)> = <S(^||)>,__sin^, 

(A8a) 

<S,(^ll)>,,_ ^{^s,-(k\\)\^y('P)\'i's.-(k\\)} = -<S(^||)>,,_cos^, 

(A8b) 

<S,(^ll)>,,_ = {"V.^Jk^Ol 2,(^) I^F^.-C^ii)) = 0. (A8c) 



Similarly, for a pure "spin-up" state I^Ps the expectation 

value of I^xi^p), ^yif), '^zif) are 



{^AkiO)s,^ =(»i',,+(^ll)|S,(^)|»I',,+(^ll)> = W||)>,_^sin^, 

(A9a) 

<S,(^ll)>,,+ =<>P,,+(^ll)|S3,(</.)|^,,+(^ll)) = -<2)(^||)),^cos^, 

(A9b) 

<S,(^ll)),,+ = <^.,-(^||)| ^xi'p) ti',s,-(k\\)) = 0. (A9c) 



(A7) 



According to Eqs. ( IA8al )-( IA9cb . one can easily find k\\ ■ 
{'L(k\\)}s^-t - 0. Therefore, for pure "spin-up" states or "spin- 
up" states, the spin orientations are in xy plane, and strictly 
perpendicular to the in-plane wave vector 
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